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Abstract  -  The  Specular-Cued  Surveillance  Web 
(SPECSweb)  multistatic  tracker  effectively  reduces  false 
track  rate  through  the  use  of  two  amplitude  thresholds. 
A  high  threshold  is  used  to  identify  the  occurrence  of 
high-strength  specular  detection  cues  for  track 
initiation.  A  low  threshold  is  used  for  selective 
extraction  of  additional  detections  for  track  extension 
forward  and  backward  in  time.  This  approach 
significantly  reduces  the  node-to-fusion-center 
communication  loading  requirements  as  well  as 
reducing  the  output  false  track  rate.  At  the  time  of  the 
specular-cued  track  initiation,  a  historical  (hack)-track 
is  reconstructed  and  output.  In  some  cases,  there  may 
still  be  false  tracks  at  the  tracker  output.  We 
demonstrate  a  post-tracking  classification  method  which 
identifies  and  removes  any  residua!  false  tracks. 
Classification  feature  information  for  track-associated 
contacts  is  passed  through  the  tracker  and  is  available  at 
its  output.  Discriminating  true  target  tracks  from  false 
ones  at  the  tracker  output  is  made  using  a  sequential 
probability  ratio  test  ( SPRT)  based  on  feature 
information  and  kinematic  track-quality.  The  method  is 
shown  to  be  effective  in  reducing  false  tracks  when 
applied  to  a  simulated  data  set. 
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1  Introduction 

A  concept  referred  to  as  the  “Specular-Cued 
Surveillance  Web  (SPECSweb)”  is  being  pursued  to 
mitigate  the  data-rate  overloading  problem  inherent  in 
multistatic  sonar  systems.  The  SPECSweb  multistatic 
tracker  does  this  through  an  implementation  of  high  and 
low  amplitude  thresholds.  The  high  threshold  is  used  to 
identify  the  occurrence  of  high-strength  specular  detection 
cues  for  track  initiation.  The  low  threshold  is  used  for  a 
selective,  proximity-based  extraction  of  additional 
detections  for  track  extension  backward  (from  the  time  of 
the  specular  cue)  in  time.  This  historical  track  data  is  then 
extended  forward  in  time  as  new  data  arrives.  This 
approach  provides  a  robust,  automated  ASW  detection 
and  tracking  method,  resulting  in  a  significant  reduction  in 
false  track  rates  and  communication  link  loading 
compared  to  conventional  multistatic  fusion  methods. 
Although  SPECSweb  has  already  been  shown  to  be  quite 
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effective  at  reducing  false  tracks  [1-3],  the  work  discussed 
in  this  paper  shows  that  yet  additional  reductions  are 
possible  through  a  post-tracking  classification  method. 

There  are  a  number  of  approaches  one  may  take  to 
perform  target  classification  together  with  sensor  fusion 
and  target  tracking,  as  depicted  in  figure  1 .  The 
traditional  approach  has  been  to  perform  classification 
prior  to  the  fusion/tracking  processes.  In  this  case,  firm 
(accept/reject)  decisions  for  detection  contacts  are  made 
prior  to  tracking.  This  acts  as  a  filter  on  the  data  being 
sent  to  the  fusion  algorithm.  It  is  usually  performed  on 
individual  contacts  from  a  single  (source-receiver- 
waveform)  “scan”,  and  classifies  each  of  them  based  on 
statistical  properties  or  physical  “features”  extracted  from 
the  data  during  upstream  signal  and  information 
processing.  The  potential  disadvantage  of  this  approach  is 
that  the  classification  process  may  not  be  accurate  enough 
or  robust,  resulting  in  unacceptable  risk  of  true  target 
dismissal.  In  this  approach,  the  success  of  tracking  will  be 
largely  dependent  on  the  effectiveness  of  the  classification 
step,  since  it  determines  which  of  the  data  are  made 
available  to  the  tracking  algorithm. 

An  alternative  approach  is  to  integrate  the 
classification  process  directly  into  the  tracking  algorithm 
itself  (sometimes  referred  to  as  “feature -aided  tracking”). 
In  this  case,  the  classification  feature  information  may  be 
used  to  more  effectively  form  tracks.  Feature  information 
may  be  used  to  aid  in  track  management  and/or  data 
association.  Likelihood  ratios  may  be  formed  on  the 
classification  feature  information,  and  incorporated  into 
track  scores,  which  can  be  used  as  the  basis  for  track 
initiation  and  termination.  The  feature  information  may 
be  used  to  improve  data  association  decisions.  The 
feature  likelihood  ratio  may  be  incorporated  into  the 
weighting  of  probabilistic  data  association  (PDA)  methods 
[4],  It  may  also  be  incorporated  into  track  scores  and 
thereby  improve  decisions  among  multiple  track 
hypotheses  (as  in  multi-frame  assignment  and  MHT 
approaches)  [5],  If  there  is  a  sufficient  ability  to  predict 
the  feature/attribute  values,  then  the  tracker  state  may  be 
augmented  to  include  the  feature  information  [4,  6]. 
These  methods  assume  that  the  feature  statistics  and  their 
distributions  (for  target  and  non-target  cases)  are  known 
well  enough  to  be  used  in  making  tracking  decisions.  If 
this  knowledge  is  not  precise  enough,  it  is  possible  that 
tracking  performance  could  be  degraded  rather  than 
enhanced. 
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A  final  approach  to  using  feature  information  with 
tracking  (and  the  one  pursued  in  this  paper),  is  to  apply 
the  classification  method  only  after  the  tracking  process 
has  been  applied.  This  allows  the  kinematic  association 
and  filtering  to  proceed  without  the  influence  of  feature 
information,  but  then  filters  the  output  based  on  the 
feature  information.  This  is  particularly  relevant  to  the 
SPECSweb  tracking  approach,  where  the  initiation 
scheme  already  significantly  reduces  the  number  of  false 
tracks  formed.  Also,  it  provides  a  backtrack  history  of 
any  of  the  target  tracks  that  it  outputs.  The  availability  of 
additional  contact  features  available  in  the  historical 
backtracks  provides  larger,  more  effective  statistical 
sample  set  for  making  certain  classification  decisions  than 
those  being  made  scan-by-scan.  Residual  false  tracks 
which  have  been  identified  may  either  be  terminated 
and/or  suppressed  from  the  operator’s  tracker  output 
display.  This  approach  will  also  be  compatible  with 
multi-sensor  classification  approaches  which  exploit  class 
features  that  depend  on  the  tracker’s  target  state  estimate. 
Future  work  will  pursue  this  extension. 

Section  2  in  this  paper  describes  the  SPECSweb 
fusion/tracking  algorithm.  Section  3  describes  the 
SPECSweb  post-tracking  classification  method.  Section  4 
shows  the  classification  method  results  on  a  simulated 
data  set.  Section  5  provides  conclusions. 


Figure  1.  Approaches  to  tracking  and  classification  (SP- 
signal  processing;  IP-  information  processing). 


2  SPECSweb  Algorithm  Description 

Detailed  descriptions  of  the  SPECSweb  multistatic 
tracking  algorithm  and  specular  cueing  approach  are 
found  in  [2-3].  A  brief  description  of  the  algorithm  is 
included  here.  In  the  SPECSweb  concept,  it  is 
envisioned  that  all  single-sensor  processing  will  be 
embedded  on  the  receiver  nodes,  and  store  all  the  local 
output  data.  Each  receiver  collects  the  contacts 
corresponding  to  one  source  transmission  as  a  single  scan 
of  “ping”  data.  Each  sonar  node  (bistatic  receiver 
processing  a  unique  source-waveform)  self-searches  each 
processed  (and  locally  stored)  scan  for  contacts  which 
exceed  a  high  SNR  threshold  setting  (HTH).  The  HTH 


identifies  very  strong  echoes,  which  likely  correspond  to 
targets  that  are  in  the  “specular  condition”.  When  in  the 
specular  geometry,  there  is  greatly  increased  target 
strength,  producing  increased  echo  energy,  as  indicated  by 
various  models  [7-8]  and  data  analyses.  The  HTH 
normally  rejects  most  (or  all)  of  the  false  alarm  clutter 
echoes,  which  have  a  lower  distribution  of  amplitudes 
than  do  specular  target  echoes  (though  not  necessarily 
lower  levels  than  non-specular  target  echoes). 

Contacts  which  cross  the  HTH  are  assumed  to  be 
“specular  cues”,  and  (only)  these  are  initially  sent  over  the 
communication  link  to  the  multistatic  fusion  center  for 
potential  track  initiation.  There  may  be  some  increased 
track  reporting  latency  using  this  approach.  Evaluation 
metrics  for  studying  the  occurrence  statistics  of  specular 
detection  in  multistatic  fields  have  been  developed  [9], 

Once  a  specular  cue  arrives  at  the  fusion  center, 
tentative  reverse-time  tracks  are  initiated.  A  cue  is 
mapped  to  an  x-y  position  in  Cartesian  coordinates,  and 
this  position  with  its  associated  error  covariance  are  then 
sent  as  a  snippet  request  to  other  nodes.  These  nodes 
calculate  the  appropriate  snippet  boundary  in  their 
respective  measurement  spaces  within  which  data 
association  would  be  possible,  according  to  a  specified 
gating  parameter.  Any  contacts  found  within  the  snippet 
gate,  and  above  the  standard  low-threshold  (LTH),  are 
then  sent  over  the  communication  link  for  further 
processing.  As  track  estimates  are  obtained,  they 
themselves  are  used  as  the  cues  for  selective  data  retrieval 
on  prior  scans  stored  on  any  of  the  nodes.  The 
retrospective  tracking  (backtracks)  continue  until  they 
meet  a  track  termination  criterion.  Recovering  track 
history  in  this  fashion  provides  valuable  contextual  and 
track  classification  information. 

The  contacts  belonging  to  the  selected  backtrack  are 
then  re-filtered  in  the  forward-time  direction,  until  the 
current  time  (of  the  initiating  specular  cue)  is  reached.  At 
this  point  the  track  classification  method  may  be  applied 
on  the  retrospective  historical  track  data.  Tracking  and 
classification  updates  continue  in  the  forward-time 
direction  updating  with  measurements  found  within  the 
retrieval  snippets  of  future  scans. 

SPECSweb  uses  logic-based  track  management  (M  of 
N  initiation  and  K  missed  detections  for  termination).  The 
tracker  is  implemented  as  an  extended  Kalman  filter 
(EKF)  with  nearest  neighbor  data  association  and  a  nearly 
constant  velocity  (NCV)  motion  model. 

3  Post-Tracking  Classification 

At  the  time  of  the  specular  initiation  cue,  a 
retrospective  backtrack  may  be  available.  The 
classification  information  within  backtracks  are 
sequentially  processed  in  a  scan/time -ordered  fashion.  The 
process  continues  as  forward-time  tracking  results  are 
newly  generated.  Classification  status  is  computed  along 
the  entire  time  duration  of  each  track.  Classification 
decisions  may  be  used  to  suppress  false  tracks  from  the 
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operator  display  (while  continuing  to  track  them),  or  to 
terminate  them. 

3.1  Sequential  Probability  Ratio  Test 

The  basis  for  the  post-tracking  classification  decision 
is  provided  by  the  sequential  probability  ratio  test  (SPRT) 
[10].  We  consider  the  two  hypotheses,  one  for  target  and 
the  other  for  non- target,  as: 


Ha  :  =>  non -target  ^ 

H  |  :  =>  target 

The  SPRT  is  given  as  the  cumulative  sum  of  the  log  of  a 
number  ( J)  of  independent  likelihood  ratios  formed  from 
the  data  representing  the  features  under  consideration  for 
classification.  At  the  nth  update  the  SPRT  is  given  by 


L(n)  =  X  In 


P{DU\HX)  P(D2j\Hx)  P(D;i  \H  j  ) 


(2) 


KP(DU\H0)  P(D2i\H0 )  P(Z>w|ff0), 

The  various  D’s  are  measurable  quantities  from  the  data. 
We  may  then  establish  two  thresholds: 


if  Hx  and  L(k )  >  A(k) 

=>  maintain  Ht 
=>set  A(k  + 1)  =  L(k) 

=>set  B(k  + 1)  =  L(k)~  A 
else  if  Ht  and  L(k)<B(k) 

=>  change  to  H0 
=>  set  A(k  + 1)  =  L(k)  +  A 
=>  set  B(k  +  l)  =  L(k) 
else  if  H0  and  L(k)>A(k) 

=>  change  to  Ht 
=>  set  A(k  +  l)  =  L(k) 

=>set  B(k  + 1)  =  L(k)~  A 
else  if  H0  and  L(k)<B(k) 

=>  maintain  H0 

=>  set  A(k  + 1)  =  L(k)  +  A 

=>  set  B(k  +  l)  =  L(k)  ^ 

end  if 

where  k  indicates  the  current  update  and  A  is  the 
difference  of  the  thresholds,  i.e., 


A  =  In 


B  =  In 


(3) 


A  =  A(X)-B(\) . 


(6) 


where  a  is  the  false  alarm  probability  (type  I  error)  and  /? 
is  the  true  target  dismissal  probability  (type  II  error),  and 
are  parameters  specified  by  the  user.  In  our  application, 
the  SPRT  is  continuously  updated  until  an  initial 
classification  is  indicated  by  a  threshold  crossing  as: 

Continue  testing  (pending)  if  B  <  L(n)  <  A 
Accept  Hx  :  =>  if  L(n)  >  A 

Accept  H0  :  =>  if  L(n )  <  B 


where  A(l)  and  B(l)  are  set  to  A  and  B  in  Eq.  (3). 

Figure  2  shows  a  comparison  of  an  example  SPRT 
output  corresponding  to  a  track  (#28  of  section  4),  which 
changes  state  as:  false  (negative  slope)  — >  true  (positive 
slope)  — *  false  (negative  slope).  The  track  is  initiated  by  a 
specular  cue  around  which  it  is  a  true  track.  Further 
backward  and  forward  in  time  the  track  updates  are 
provided  by  false  contacts  and  the  track  becomes  false. 
The  adjustable  threshold  method  makes  an  earlier  and 
better  classification  of  the  target  section  than  the  fixed 
method  (scans  2301-3250  vs.  scan  2900-3250). 


We  also  wish  to  consider  the  possibility  that  some 
tracks  may  change  status,  i.e.,  change  from  being  a  true 
target  track  to  becoming  a  false  track,  or  vice  versa.  This 
can  occur  in  tracking  targets  amongst  dense  clutter,  where 
a  true  track  starts  associating  nearby  false  contacts  and 
then  progressively  “wanders”  off.  In  order  to 
accommodate  this  possibility  we  have  chosen  to 
implement  adjustable  thresholds  after  the  initial 
target/non-target  determination  has  been  made.  This  is 
done  by  resetting  the  thresholds  by  raising/lowering  them 
according  the  new  maximal/minimal  values  of  L.  The 
following  adjustment  and  testing  scheme  is  used: 


Figure  2.  SPRT  (black)  with  fixed  thresholds  (top)  and 
adjustable  thresholds  (bottom);  target  threshold  (red),  non¬ 
target  threshold  (blue). 


3.2  Kinematic  Track  Quality 

Many  tracking  algorithms  utilize  an  SPRT  for  track 
initiation  which  includes  both  a  kinematic  and  a  signal- 
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related  contribution  [11].  Here  we  use  only  the  kinematic 
portion  of  this  approach,  and  it  will  be  applied  to  the 
tracker  output  rather  than  the  tracker  input.  The  signal- 
related  contribution  depends  on  a  priori  knowledge  of  PD 
and  PFA,  which  in  practice  may  be  difficult  to  know  well 
enough.  While  PFa  may  be  ascertained  by  analysis  of  in- 
situ  data,  PD  for  multistatic  systems  is  complex  to  estimate, 
being  dependent  on  a  complex  space/time  evolving 
acoustic  environment,  as  well  as  the  source-receiver 
positions,  target  aspect  angle,  etc.  These  dependencies 
may  be  estimated  using  sophisticated  sonar  equation 
models;  however  they  are  often  still  not  of  sufficient 
fidelity  to  adequately  predict  precise  performance.  An 
approach  to  estimating  PD  from  the  arriving  data  has  been 
proposed  [12].  Nevertheless,  given  the  difficulties  of 
predicting  PD,  we  choose  to  focus  on  only  the  kinematic 
contribution  in  this  work. 

The  assumption  for  the  kinematic  likelihood  ratio  is  that 
target-originated  contacts  will  be  distributed  according  to 
Gaussian  statistics  with  the  mean  at  the  state  estimate, 
whereas  false  alarm  contacts  are  assumed  to  be  uniformly 
distributed  over  the  association  gate.  The  likelihood  ratio 
for  the  kinematic  portion  of  the  SPRT  is  given  by  [  1 1  ] 


LRk  = 


(7) 


where  the  volume  of  the  ellipsoidal  association  gate  is 
given  by 


Vc  = 


7TM'2-GM/1-J\S\ 


(8) 


r  is  the  gamma  function,  M  is  the  dimensionality  of  the 
measurement  and  G  is  the  value  corresponding  to  a  gate  of 
a  specified  association  probability  (depends  also  on  M).  S 
is  the  Kalman  residual  covariance  matrix,  given 
by  S  =  HPHt  +  R  ,  where  H  is  the  linearized  measurement 
(Jacobian)  matrix,  R  is  the  measurement  covariance 
matrix,  and  P  is  the  predicted  state  covariance  matrix,  d2 
is  the  norm  of  the  residual  vector  given  by 


when  there  is  no  contact  associated  with  the  track  update 
(i.e.,  a  track  “coast”  condition). 

3.3  Classification  Feature  Score 

We  assume  that  per-contact  classification  feature  scores 
have  been  computed  in  the  processing  upstream  of  the 
tracker.  These  may  be,  for  example  echo  SNR-based  or 
related  echo  energy  size/shape.  Here,  we  don’t  address 
any  specific  or  particular  feature,  but  consider  the  case  of  a 
“generic”  feature,  for  which  a  score  is  computable.  When 
a  contact  is  assigned  to  a  track  update,  information  about 
that  contact,  including  the  feature  score  is  saved  and 
output  along  with  the  updated  state  estimate. 

We  assume  that  we  have  some  knowledge  of  the 
underlying  target  and  non-target  distributions  for  the 
generic  feature  score  that  is  computed.  If  these 
distributions  are  not  known,  target  and  non-target  training 
data  may  be  used  to  obtain  estimates  of  their  probability 
density  functions  using  kernel  density  estimation  methods. 
Given  the  target/non-target  distributions,  and  the  feature 
score,  we  can  form  the  likelihood  ratio  for  the  feature  as 

LRf 

which  is  a  ratio  of  the  corresponding  probability  densities, 
given  the  value  of  the  feature  score  DF. 

The  likelihood  calculations  are  made  after  each  update 
of  the  Kalman  filter  if  a  contact  updates  the  track  estimate. 
If  the  track  is  in  a  “coast”  condition,  then  the  likelihood 
ratio  is  not  formed  because  there  is  no  associated  feature 
information.  The  SPRT  is  updated  sequentially  as  new 
data  becomes  available. 

3.4  Combined  Classification 

The  post-track  classification  may  utilize  both  of  the 
likelihood  ratios  (kinematic  and  feature)  just  described.  In 
this  case  the  SPRT  becomes: 

L{n)  =  f\n[LRK(i)LRF(i)}  (12) 

i=l 


p(PF  Igi) 
P(Df\H0) 


d1  =  \v{k)-h(x(k\k-\)Y  S-'\v(k)-h(x(k\k-\)\  (9) 


4  Classification  Results 


where  v  is  the  measurement  data,  x(k\k  - 1)  is  the 

projected  target  state,  and  h  is  the  non-linear  measurement 
function.  Simplifying  (5)  yields 


_  g(_rf2/2)  ■  (G/2)m/2 

K~  r((M/2)  +  l) 


(10) 


The  calculations  to  obtain  d2  are  made  during  iterations 
of  the  Kalman  filter  update,  and  are  saved  and  output 
along  with  the  tracker’s  updated  target  state  estimates. 
The  likelihood  ratios  are  formed  and  the  SPRT  is 
calculated  at  each  iteration  of  the  Kalman  update,  except 


Multistatic  data  was  simulated  using  the  Passive- 
Active  Contact  Simulation  Tool  (PACsim)  [13],  Figure  3 
shows  the  scenario,  with  5  fixed  sources  (red)  and  16 
fixed  receivers  (green)  distributed  over  a  surveillance 
area.  Pings  are  emitted  every  minute,  cycling  amongst  the 
sources,  and  alternating  between  FM  and  CW  waveforms. 
The  scenario  lasts  about  6.8  hours,  with  409  pings  and 
6,541  measurement  scans  (16  scans  per  ping).  Eight 
targets  are  simulated  as  described  in  Table  1  and  as 
depicted  in  figure  3. 
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Figure  3.  Simulation  Scenario  (distances  in  10  km); 
sources  (red),  receivers  (green),  targets  (black). 


Table  1.  Simulated  Targets 


# 

Description 

1 

6  course  changes  of  ±45°;  (loose  zig-zag),  5  kts 

2 

Fixed  target  to  the  north 

3 

6  course  changes  of  1 80°;  back  and  forth  5  kts 

4 

Fixed  target  to  the  south 

5 

6  course  changes  of  90°  ;  (box),  7.5  kts 

6 

Fixed  target 

7 

6  course  changes  of  135°;  (tight  zig-zag),  5  kts 

8 

Fixed  target  to  the  east 

The  target  SNRs  are  computed  by  solving  the  sonar 
equation  for  each  so  u  re  c  -  rcc  c  i  vc  -  target-  wa  vc  fo  r  m 
combination  in  a  representative  acoustic  environment  with 
generic  airborne  multistatic  system  components.  The 
simulation  includes  modeling  of  transmission  loss, 
reverberation,  aspect-dependent  target  strength,  noise 
background,  and  system  gains.  Measurements  of  contact 
time  delay  (r),  bearing  ((9),  and  Doppler  (v)  were 
simulated  using  assumed  errors  of  o  =0.1  sec  (FM), 
a  =0.4  sec  (CW),  Og=4°  and  a, ,=0.4  m/s.  Non-target 
(clutter)  contacts  were  distributed  uniformly  in 
measurement  space  and  with  SNRs  following  an 
exponential  distribution  above  the  detection  threshold. 
The  number  of  false  contacts  per  scan  is  Poisson 
distributed  (2=30)  per  scan.  A  generic  classification 
feature  score  was  simulated  for  each  contact  with  random 
draws  from  the  following  normal  distributions: 

•  Non-Target:  N(jUc  =  0,  <7C  =  10) 

•  Target  (easel):  N(jUn  =  6, cr,  =  10) 

•  Target  (case  2):  N(jUT2  =  10,  a,  =10) 

Figure  4  shows  a  plot  of  these  distributions  and  their 

associated  ROC  curves  vs.  the  chance  line  (indicated  in 
green).  We  see  that  target  case  2  has  more  separation 
from  the  clutter  distribution  than  case  1,  and  therefore, 
greater  classification  potential.  Figure  5  shows  the 
cumulative  output  of  the  SPECSweb  tracker  with  a  good 
choice  of  parameter  settings.  We  see  that  the  result  for 
this  case  is  excellent:  good  tracking  with  no  false  tracks. 
The  average  input  and  output  PD  were  0.3  and  0.9, 
respectively.  The  target  tracks  average  purity  (target-to- 


total  contact  updates)  was  91%  and  the  percentage  of 
“coasts”  (scans  without  a  contact  update)  was  92%.  For 
this  scenario,  performance  is  so  good  that  it  does  not 
require  post-track  classification. 


Figure  4.  Generic  classification  feature  score  (Target 
easel-blue,  Target  case2-red,  Clutter-black);  distributions 
(left),  ROC  Curves  (right). 
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Figure  5.  Optimized  SPECSweb  tracking  output  without 
classification;  target  true  trajectories  (yellow),  tracks 
(blue),  sources  (red),  receivers  (green). 


(meters) 


Unlike  the  performance  just  seen,  more  challenging 
scenarios  and  acoustic  environments  may  yield  a  number 
of  false  output  tracks,  even  with  the  SPECSweb 
algorithm.  In  order  to  create  such  a  test  data  set  for  our 
classification  analysis,  we  now  select  sub-optimum 
tracking  parameters  to  induce  an  increase  in  the  false  track 
rate.  This  was  done  by  changing  the  track  management 
logic  to  make  it  easier  to  initiate  and  harder  to  terminate 
tracks  (M/N/K  =  6/1 6/50  vs.  M/N/K=  6/81/150).  Figure  6 
shows  the  resulting  tracker  performance.  We  see  a 
cluttered  output  display,  composed  of  65  tracks,  some  of 
them  extending  beyond  the  surveillance  box  boundary. 
This  test  data  will  serve  as  a  baseline  case  (without 
classification)  and  as  the  input  to  the  post-tracker 
classifier.  Some  tracks  are  true  target  tracks;  others  are 
initiated  by  true  target  contacts  and  continue  on 
(wandering)  with  false  alarms,  or,  they  are  wholly  and 
completely  false  tracks.  The  mapping  of  targets  to  output 
track  number  is  given  as:  TGT1  (10),  TGT2  (2,58),  TGT3 
(4,7),  TGT4,  (1,31,59,64),  TGT5  (3),  TGT6  (8,28,38), 
TGT7  (5),  TGT8  (11,39). 


The  classification  method  is  now  applied  to  the  tracker 
output  (of  figure  6).  The  first  case  we  consider  is  to  use 
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only  the  kinematic  information,  as  described  in  section 
3.2,  with  decision  thresholds  set  corresponding  to 
a=P=0.001.  Figure  7  shows  the  entire  time-extent  of  each 
of  the  65  input  tracks  (referenced  by  scan  index),  from  the 
earliest  point  of  the  backtrack  to  the  last  update  of  the 
forward  track  termination.  The  time  of  the  track-initiating 
specular  cue  is  shown  overlaid  by  the  magenta  marker. 
The  tracks  are  color-coded  by  the  classification  status 
obtained:  target  (red),  non-target  (green),  and  decision 
pending  more  information  (blue).  We  observe  that  many 
of  the  tracks  are  classified  as  false,  and  that  the 
classification  decisions  are  usually  made  prior  to  the 
occurrences  of  the  specular  cues  (tracker’s  initial  output). 
Most  of  the  target  tracks  are  detected  early  and  have 
continuous  hold,  except  in  a  few  instances  when  the 
classification  status  changes  from  true-to-false  (or  vice 
versa).  This  is  due  to  tracks  which  “wander”,  being  fed  by 
false  alarms  and  which  most  often  occur  for  the  fixed 
(non-mobile)  targets.  The  true  (red)  and  pending  (blue) 
tracks  are  shown  plotted  in  figure  8,  where  we  see  good 
correspondence  between  the  true  tracks  and  the  eight 
known  target  trajectories,  and  a  significant  reduction  in  the 
false  (clutter)  tracks  observed.  The  classifier  in  this  case 
(figure  8)  shows  good  performance  compared  to  the  input 
case  (figure  6). 


x  104 

Figure  6.  Data  set  with  induced  false  tracks  for 
classification  evaluation;  65  output  tracks  (blue),  target 
true  trajectories  (yellow). 

Next  we  show  examples  only  using  the  generic  feature 
information,  as  described  in  section  3.3,  with  decision 
thresholds  nominally  set  corresponding  to  <*=/?=  0.001. 
Figure  9  shows  the  output  of  the  classifier  for  Case  1 
(heavily  overlapping  target/non-target  distributions).  Here 
it  is  observed  that  classification  decisions  take  a  longer 
time  to  be  reached,  as  indicated  by  the  increases  in  the 
amount  of  pending  status  (blue).  Although  the 
classification  decisions  for  targets  are  quite  good,  many  of 
the  false  tracks  are  not  classified  prior  to  the  specular  cue 
initiation.  Figure  10  shows  the  track  display  for  this  case, 
which  shows  good  tracking  on  the  targets,  but  with  a  few 
falsely  classified  true-tracks  appearing.  There  is  a  large 


number  of  pending  tracks  shown  which  confuse  the 
output.  One  could  choose  to  display  only  confirmed  true 
target  tracks,  if  desired,  which  change  the  impression 
significantly.  Nonetheless,  the  results  are  still  much  better 
than  the  input  case  (figure  6). 

Figure  1 1  shows  the  effect  of  relaxing  the  false  alarm  and 
false  dismissal  parameters  (changing  a=fl=  0.001  to 
a=P=  0.01).  Classification  decisions  are  reached  much 
faster  (fewer  observable  pending  tracks),  as  expected, 
though  perhaps  with  worsening  false  track  rate. 
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Figure  7.  Classification  vs.  track,  for  kinematic 
information  only;  pending  (blue),  non-target  (green), 
target  (red);  specular  track  initiation  (magenta). 
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Figure  8.  Output  for  kinematic  information  only;  target 
trajectories  (yellow),  classified  target  (red);  pending 
classification  (blue). 
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Figure  9.  Classification  vs.  track,  for  feature  information 
only  (case  1);  colors  as  in  figure  7. 
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Figure  10.  Classification  results  for  feature  information 
only,  case  1,  a=/2= 0.001. 
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Figure  11.  Output  for  feature  information  only,  case  1, 
a=P=0.01. 


Figures  12  and  13  correspond  to  the  Case  2,  where  the 
feature  score  is  more  separable  (and  with  original  settings 
of  a=/?=0.00 1 ).  Here  we  obtain  even  better  performance, 
though  it  is  not  as  good  as  the  results  obtained  with 
kinematic  information  only.  Here,  the  total  percentage  of 
pending  track-time  (of  total  track-time)  is  reduced  to  17% 
from  42%  (see  in  figures  9). 

Next  we  consider  the  effect  of  combining  the 
kinematic  and  feature,  as  described  in  section  3.4,  with 
decision  thresholds  again  corresponding  to  a=/3== 0.001. 
The  results  are  shown  in  figures  14.  In  this  case  we  obtain 
the  best  results,  with  the  cleanest,  de-cluttered  output 
display.  All  eight  targets  are  tracked  well.  When  true 
tracks  wander  away  from  real  targets,  they  are  classified 
false  and  suppressed.  Classification  decisions  are  made 
earlier,  and  changes  in  state  are  more  robust  to  the  truth. 
This  example  shows  the  power  of  combining  all  available 
information  into  the  classifier  likelihood  ratio  formulation. 


Scan  #) 

Figure  12.  Classification  vs.  track,  for  feature  information 
only  (case  2);  colors  as  in  figure  7. 


Figure  13.  Output  for  feature  information  only,  case  2, 
«=/?== 0.001. 
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Figure  14.  Output  for  combined  kinematic  and  feature 
information,  case  2,  a=f}= 0.001. 

5  Conclusions 

A  post-tracking  classifier  method  compatible  with 
SPECSweb  has  described,  and  has  been  shown  to  be 
effective  in  reducing  false  tracks  in  the  output.  So  far,  the 
formulation  includes  kinematic  track-quality  measures 
and  per-contact  feature  scores.  Their  likelihood  ratios 
have  been  used  in  an  adjusting  sequential  probability  ratio 
test  to  arrive  at  classification  decisions.  Often,  these 
decisions  are  made  during  the  backtrack  section  of  the 
total  track  (before  the  specular  cue)  and  the  non-target 
and/or  pending  tracks.  Table  2  summarizes  the 
performance  of  the  cases  evaluated.  Future  work  will 
investigate  the  incorporation  of  specific  target-aspect  (and 
therefore  tracker  output  state  estimate)  dependent  field 
features. 


Table  2.  Out] 

out  class  percentages 

# 

% 

True 

% 

False 

% 

Pend 

Kinematic 

21 

66 

13 

Feature  (case  1) 

25 

33 

42 

Feature  (case  2) 

27 

56 

17 

Both  (case  1 ) 

22 

66 

11 

Both  (case  2) 

23 

67 

9 
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